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ABSTRACT 

Sixty percent of gamma-ray bursts (GRBs) reveal strong Mgll absorbing systems, which is a 
factor of ^2 times the rate seen along lines-of-sight to quasars. Previous studies argue that the 
discrepancy in the strong Mgll covering factor is most likely to be the result of either quasars 
being obscured due to dust, or the consequence of many GRBs being strongly gravitationally 
lensed. We analyze observations of quasars that show strong foreground Mgll absorption. We 
find that GRB lines of sight pass closer to bright galaxies than would be expected for random lines 
of sight within the impact parameter expected for strong Mgll absorption. While this cannot be 
explained by obscuration in the GRB sample, it is a natural consequence of gravitational lensing. 
Upon examining the particular configurations of galaxies near a sample of GRBs with strong 
Mgll absorption, we find several intriguing lensing candidates. Our results suggest that lensing 
provides a viable contribution to the observed enhancement of strong Mgll absorption along 
lines of sight to GRBs, and we outline the future observations required to test this hypothesis 
conclusively. 

Subject headings: cosmology: Gamma ray bursts, gravitational lensing 



1. Introduction 



Porciani et al. (2007) considered dust obscura- 



Prochter et al.l ( 2006 ) pointed out that gamma- 



ray bursts (GRBs) show approximately four times 
as many Mgll absorbing lines in their spectra as 
quasars, although the current best estimate of the 
excess is now a factor of ~2 ( [VerganT "eral][2009l 
hereafter V09). Strong Mgll absorbers (equiva- 
lent width, EW, greater than lA) are found to- 
wards ^60% of GRBs with foUowed-up optical af- 
terglows. Both being high redshift beacons, GRBs 
and quasars might be expected to have similar 
lines of sight through the cosmos, and explain- 
ing the preponderance of Mgll absorbers towards 
GRBs has proven a challenge 
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tion, beam size differences, the intrinsic proper- 
ties of GRBs, and gravitational lensing as possible 
causes. Beam size differences were found to be ir- 
relevant, with simulations predicting that the ab- 
sorbing systems are significantly larger than either 
the GRB afterglow ring or the quasar accretion 
disk. If the additional strong Mgll absorbers are 
physically associated with the GRBs, one might 
expect the absorber properties to be distinct from 
the systems towards quasars. However, Cucchiara 



et al. (20091 found no difference between the two 
populations of absorbers. Ultimately, this leaves 
dust obscuration of quasars and gravitational lens- 
ing of GRBs as the two most plausible explana- 
tions. 

If strong Mgll absorbers are dusty, so that the 
discrepancy between GRBs and quasars is due to 
quasars being preferentially lost from the observed 



for either weak Mgll systems ||Tejos et al.|2009| l or for CIV 
absorption systems | |Tejos et al.|2007^ . 
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samples (GRB's being initially selected in gamma- 



rays are not subject to this extinction), Porciani 



et al. ( 2007 ) determined that the number of Mgll 



systems towards quasars would have to have been 
underestimated by a factor of 1.3-2, which they 



found to be unlikely (but see Budzynski & Hewett 
2011). Studies of quasars have shown that the 



metal-enriched gas responsible for Mgll absorp- 
tion is broadly associated with foreground galax- 
ies (e.g. Kacprzak et al. 2007). Thus, if dust 



obscures quasars with strong Mgll, providing the 
origin of the GRB-quasar absorption system dis- 
crepancy, then the ~60% of GRBs with strong 
Mgll absorbers should represent random lines of 
sight within 60% of the sky nearest to foreground 
galaxies. However, if lensing is responsible for the 
discrepancy of Mgll absorption between quasars 
and GRBs, then there should be an excess of 
GRBs at small separations from foreground galax- 
ies. Therefore, the distribution of GRB-galaxy 
separations provides a test to distinguish between 
the two hypotheses. 

Finding multiple images of the same source is 
the calling card of strong gravitational lensing. 



Wyithe et al.||2011| (Wll) explored the possibility 
that gravitational lensing combined with a multi- 



band magnification bias (described in Wyithe 
et al.|2003 ), could lead to the large number of high 



equivalent width Mgll absorbers for GRBs. Wll 
showed, if the gamma ray and optical luminosities 
of GRBs are uncorrelated, and if the luminosity 
functions have a cumulative slope with power law 
index > 3.5, then 10%-60% of the GRBs with 
afterglow follow-ups should have been multiply- 
imaged. 

The V09 sample was chosen solely on the basis 
of the optical afterglow's brightness and reveals 15 
strong Mgll absorbers out of the 26 GRBs sample. 
As Wll do not predict lensed GRBs to be signif- 
icantly brighter overall, we treat the V09 sample 
as an unbiased group of GRBs and therefore ex- 
pect 10%-60% of the 26 GRBs should have been 
multiply-imaged. Wll also noted that, based on 
the sky coverage of Swift, the probability of the 
satellite detecting three separate doubly-imaged 
GRBs is only 3%. The probability that none of the 
V09 sample was observed to be doubly-imaged, 
even if lensing does occur, was found to be ^50%. 

If GRBs are frequently strongly lensed, they 
open up a potentially powerful probe of galax- 



ies and the cosmos. Time delay and magnifica- 
tion measurements, as done with quasar lenses, 
become possible with exquisite precision, allowing 
accurate mass models and distance measurements 
to be made to an ever increasing sample of objects. 

In this paper, we investigate the possibility that 
gravitational lensing is affecting a significant frac- 
tion of GRBs with optical counterpart. For con- 
sistency with the Wll analysis, we study archival 
Hubble Space Telescope (HST) observations of the 
galaxies near the lines of sight towards each of 
the GRBs in the V09 sample, and estimate the 
probability that the GRB was strongly lensed. 
In section [2j we describe the data used and our 
data reduction procedures. In section [3j we study 
the probability of finding nearby galaxies towards 
GRBs with strong Mgll absorbers vs. a random 
line of sight. In section |4j we outline the lensing 
analysis methods and properties of the model em- 
ployed, and present our results for the individual 
GRBs. We discuss our conclusions in section |5] 
Throughout this paper we assume a cosmology of 
nM = 0.27, r^A = 0.73 and -ffo=72 km/s/Mpc. 

2. Data 

In order to consistently test the Wll hypoth- 
esis, we analyze the same 26-GRB sample from 
V09, from which their model is based. As we are 
searching for galaxies with small impact parame- 
ters, we require high spatial resolution, which is 
best obtained with HST. In the few cases where 
HST data are not available, we use data from 8-m 
class telescopes. Due to the relatively poor an- 
gular resolution of the VLT and Gemini imaging 
data, a null detection in these cases does not rule 
out lensing. These data are taken as a part of our 
statistics only when a possible lensing galaxy is 
observed. 

Calibrated HST images available from the 
archive, taken with the STIS and AGS instru- 
ments, were combined in python by way of the 
STSB AS /MuUidrizzle package^ VLT ISAAC, 
FORSl and F0RS2 images were similarly reduced 
using the relevant pipeline within the (7as(?anc[^ 
software package. Version 1.9 of the Gemini 



^http: / /www.stsci.odu/hst/HSToverview/documents/multidrizzle 
•^http: / /www. cso.org/sci/software/gasgano/ 
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IRAF packag^ was employed in reducing data 
from Gemini's GMOS instrument (Hook et al 



2004). After stacking all images from each filter 



and constructing deep images from our reduced 
data, nearby galaxies (< 5") were located using 
SExtractor (iBertin & Arnoutsl 119961). The im- 



pact parameter between these galaxies and their 
partner GRBs was subsequently measured. 

In total, we have deep images of 11 objects with 
strong Mgll absorbers from the sample of V09 
(representing ~75% of the 15 GRBs with strong 
Mgll absorbers). Since the HST observations were 
not made on the basis of strong Mgll absorption, 
we believe that this is an unbiased sample for this 
analysis. 

3. Galaxy Proximity to GRB Lines-of- 
Sight 



Ghen et al. (2009) noted the consistent pres- 



ence of galaxies at small angular separations from 
a sample of four GRBs showing strong Mgll ab- 
sorbers, and the absence of such galaxies near 
three GRBs without them. Such a correspon- 
dence is not surprising since the Mgll absorption 
is thought to arise from material associated with 
galaxies. Our analysis takes this approach fur- 
ther by roughly tripling the sample size, and by 
comparing Mgll lines of sight to the more robust 
baseline of a distribution constructed from random 
lines of sight. 

According to V09, ~60% of all GRBs have a 
strong Mgll absorbing system (> lA). Using the 
galaxy catalog of Coe et al. ( 2006[) for the Hu bble 
Ultra Deep Field (UDF; [Beckwith et al.|2006| ) , we 
measure the radius around every galaxy brighter 
than a given magnitude which, together, would 
give 60% sky coverage using a top hat function; we 
denote this radius as OmqII- This top-hat model 
represents the most extreme concentration pos- 
sible for Mgll absorption along random lines of 
sight, and thus is the most stringent comparison 
for the observed impact parameter distribution. 
It may not, however, be an unrealistic model, as 
Ghen et aT] ( |2010 1 used quasar sight-lines to mea- 
sure a high Mgll covering fraction within a certain 
(luminosity-dependent) galaxy radius, and found 
that the covering fraction fell sharply at larger 

*http:/ /www. gemini.edu/sciops/ 



radii. 

If GRBs are randomly distributed, one would 
expect the 60% of GRBs which have Mgll absorp- 
tion to have a random distribution of separations 
from galaxies within the radius OmqII- However, if 
the cause for the discrepancy between GRBs and 
quasars is due to gravitational lensing, we would 
expect the distribution of galaxies around GRBs 
with strong Mgll lines to be more concentrated 
towards small radii. In the framework of our UDF 
analysis, the value of 0Mgii in any filter is a func- 
tion of only the limiting magnitude used to con- 
strain the galaxy catalog. 

3.1. Distribution of separations to nearby 
galaxies 

We use a Monte Carlo approach to generate 
random lines of sight through the UDF, and mea- 
sure the distances to the nearby galaxies in the 
Goe et al. catalog. We adopt a limiting magnitude 
for the catalog galaxies of 27.5 mag (AB) in each 
filter, since this corresponds to the faintest nearby 
galaxy in the GRB images. Running 10,000 real- 
izations in each of the four UDF filters, we then 
weight the relative contributions from each filter 
to match the filter distribution for the GRB im- 
ages (using only the filter that gives the deepest 
stacked image for each GRB in the analysis, and 
taking the UDF filter most similar to that used 
with the GRB). 

Our statistical analysis does not distinguish be- 
tween galaxies that lie in front of the GRBs and 
those behind. As long as the galaxy number 
counts continue to rise at the faint end (a con- 
dition satisfied for F775W=27.5 mag [Fig. 29 of 
Goe et al. 2006], as well as for the other UDF 
filters), the consequence of our blindness to the 
galaxy redshifts is to add spurious objects to our 
analysis for both the GRBs and the random lines 
of sight. This will weaken the signature of lensing 
by galaxies at small impact parameters by mixing 
in unrelated background objects and reducing the 
difference in the cumulative radial distributions. 
If Mgll absorption is related to galaxies that are 
typically fainter than any of our images, then we 
would expect the spatial distribution of galaxies to 
be the same for the GRBs and the random lines of 



sight. Bornancini et al. (2004) and Gampisi et al 



(2009) found that long-GRB host galaxies popu- 



late lower density regions than average, and we 
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therefore do not expect a bias due to galaxies at 
the redshift of the GRB. 

A Kolmogorov-Smirnov (K-S) two-sample test 
rejects the hypothesis of the GRB data being 
drawn from the random distribution, but only at 
> 90% confidence when using a limiting magni- 
tude of 27.5; because we would not expect the 
entire population to be lensed, but rather a few 
very unusual cases, it is not surprising that this 
treatment is inconclusive. 

3.2. Probabilities of nearest bright galeLxy 

Assuming that lensing is caused by the nearest 
galaxy, we perform a test to compare the prob- 
ability of the 11 different alignments between the 
GRB and its nearest galaxy. This time, we find the 
probability of each event happening when compar- 
ing to a line of sight near a galaxy with the same 
apparent brightness as our GRB host. While, as 
mentioned above, GRB host galaxies are not likely 
to reside within a cluster, this is done to account 
for any improbable clustering, and falsely identi- 
fying galaxies that are associated with the GRB 
redshift itself. 

After randomly choosing a galaxy from the 
UDF with similar apparent magnitude as the ac- 
tual GRB host (±0.5 mags), we measure the mag- 
nitude of and distance from it to the nearest 
galaxy (up to a specific limiting magnitude). For 
GRBs without known host galaxies we use a ran- 
dom line of sight. The probability of finding a 
closer and/or nearer galaxy is calculated and the 
probabilities of the 11 alignments are multiplied 
together to form one probability for the ensam- 
ble. It is important to note that only alignments 
with a galaxy within OmqII are counted (as we as- 
sume those galaxies are the cause for the strong 
Mgll absorption). In figure [l] the histogram of 
the multiplied probabilities is plotted (for limiting 
magnitude of 27.5 in V band). The red line rep- 
resents the probability of the GRBs alignments. 
We find the probability of having a set of align- 
ments as close and as bright as the observed sam- 
ple of 11 GRBs is only 0.2% (upper panel). To 
verify that random lines of sight are representa- 
tive of the GRB population without strong Mgll 
absorption, we compare the random lines of sight 
result to that of 9 GRBs without Mgll absorp- 
tion which have publicly available HST observa- 
tions (060526, 061201, 070721B, 060313, 050730, 



050419A, 050525, 060614 and 060729). The prob- 
ability of finding the alignments from these 9 
GRBs to their nearest galaxy is 78%, confirming 
that our MC routine is adequate. 

As we remove our probable candidates for 
strong lensing and possibly magnified GRBs from 
the analysis (see Q, we expect the alignments to 
be similar to our expectation from a random line 
of sight. Therefore, we test the outcome when 
ignoring the most likely candidates GRB020405, 
GRB030429 and GRB991216. The same simula- 
tion now finds the alignment probability of the 8 
remaining GRBs to be 15% (lower panel). As the 
distribution appears to have a log-normal shape, 
P=0.15 is just within 1 standard deviation of the 
mean (P=0.5) for a random line of sight. Fol- 
lowing the statistical approach of V09, we find 
that the absorber density per unit redshift (^) 
changes from 0.74 ± 0.20 to 0.58 ± 0.20 when re- 
moving GRBs 020405, 030429 and 010222. When 
comparing to the QSO result of 0.278 ± 0.010 



(Nestor et al. 2005), the significance of the Mgll 
discrepancy is reduced by one sigma. However, 
GRBs are not found to completely agree with the 
absorbers distributions in QSO, implying another 
process (e.g. dust obscuration) must be involved. 

In order to test the effect of the limiting mag- 
nitude, we repeat the analysis for limiting magni- 
tudes in the range 27.0 - 28.5. The brighter limit 
is chosen to match the faintest nearby galaxy ob- 
served for the GRBs, and the fainter is chosen as 
a highly conservative limiting magnitude for our 
images (we do not have such deep images, how- 
ever a fainter limiting magnitude will make Omgii 
smaller, which makes it less likely to find a galaxy 
as bright or as close as the nearest galaxies to the 
GRBs). The results for the probability of the 11 
GRBs happening range from 0.1-0.4% (and for ig- 
noring the likely cases range from 12-16%), sug- 
gesting that our analysis is not highly sensitive to 
the limiting magnitude. This confirms that GRBs 
have brighter, closer galaxies than expected if the 
60% of GRBs with Mgll represent random lines 
of sight within the 60% of the sky closest to fore- 
ground galaxies. This is as expected under the 
gravitational lensing hypothesis. 

To verify our results are not highly dependent 
on cosmic variance, we perform another test. Tak- 
ing the extreme case (galaxies with ~ lO^^M©) 
for the UDF from Moster et al. (2011 1, we assume 
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30% variance within the whole UDF field. In a 
Monte-Carlo simulation we increase the number of 
galaxies in the field drawing from a gaussian with 
a mean taken from the size of the UDF catalog 
and 30% variance. Each extra galaxy is randomly 
placed in the field and is given a magnitude which 
is randomly drawn from the UDF magnitude dis- 
tribution. Repeating the process above we find 
the probability of the 11 GRBs having their align- 
ments with their nearest galaxy to increase to only 
3%. We therefore conclude that our results are not 
particularly sensitive to cosmic variance. 




0.4 0.6 
Relative Probability 



Fig. 1. — Monte-Carlo analysis describing the 
probability of having the alignments between the 
GRBs and their nearest galaxy. The cumulative 
fraction is found using a Monte-Carlo approach, 
and choosing the nearest 11 galaxies towards ran- 
dom galaxies having the magnitudes of the GRB 
hosts (±0.5). Where a GRB host is not observed, 
a random line of sight is chosen. The upper panel 
is for the 11 GRBs, showing a probability of 0.2%. 
The lower panel shows the result after remov- 
ing the 3 most likely candidates: GRB020405, 
GRB030429 and GRB991216. The probability of 
having the alignment of the 8 GRBs left is 15%. 



4. Gravitational Lensing Models of indi- 
vidual GRBs 

We showed above that GRBs with strong Mgll 
absorbers include outliers with improbable small 
separations from foreground galaxies. We there- 
fore conclude that lensing is a likely cause of the 
Mgll absorber excess since it is consistent with 
such an alignment, whereas obscuration of quasars 



is not. Armed with this motivation, we investigate 
the individual cases to look for potential lenses. 
In this section we proceed under the assumption 
that the galaxy proximity is related to gravita- 
tional lensing of the GRBs and conduct detailed 
studies of the lensing likelihood for each GRB in 
the V09 sample. 

It is important to note that like our statistical 
analysis, the examination of lensing for individual 
GRBs may be confused by the presence of galax- 
ies lying behind the GRB in question. While we 
can use the UDF redshift catalog of Coe et al. 
(2006) to constrain the likely distances to some of 
our galaxies based on their SEDs, determination 
of the actual redshifts can only be accomplished 
by future observations. 

We determine the properties of all potential 
lensing galaxies using the methods of Tucker et 
al. (2012a, b in prep). In summary, galaxy spec- 
tral energy distribution (SED) templates are cre- 
ated from the UV to Infrared using 11 galaxies 
that cover a complete range of morphological types 



(Calzetti et al. 


|1994 


1996 




Kinney et al. 


1996 


Mannucci et al.| 2001 


. We reproduce the range 



of observed galaxies via linear combinations of all 
SEDs, and find all acceptable fits to the observed 
photometry using a statistic. From the range 
of acceptable SEDs, we calculate rest-frame mag- 
nitudes and the x^-weighted galaxy classification 
when data in several bands is available. 

In the few cases where we identified the galaxy 
as a spiral and could estimate M^, we use the 



Tully-Fisher relation (Miller et al. 2011) to ap- 



proximate the rotational velocity of stars in the 
lensing galaxy. In order to convert from the rota- 
tion al velocity to velocity disp ersion, we divide by 
\/2 ( Binney fc Tremaine[l987 |. Incorporating the 
angular impact parameter, the estimated velocity 
dispersion, and the known absorber redshifts into 
a singular isothermal sphere (SIS) model, we use 
the observational and model uncertainties to cal- 
culate the probability of magnification and of the 
lens producing two images of the GRB. 

We use the GRA VLENS/LENSMODEL soft- 
ware (Keeton 20011 to model the complex sys- 



tems (i.e., those with multiple lensing galaxies). 
The software allows one to specify the locations 
of the source (GRB), lenses, and images of the 
source (including time delay, when applicable), to 
constrain possible mass profiles of the lenses. The 
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output includes the image magnifications and time 
delays (when not provided) , and the critical curve 
of the lenses in the image plane. In the few cases 
where more than one galaxy is observed near the 
GRB and where none of the galaxies is obviously 
associated with the absorbing system (via a spec- 
troscopic rcdshift for the galaxy or extreme close 
proximity to the GRB), we model them as a group 
where each galaxy is associated with a Singular 
Isothermal Ellipse (SIE) with shear. Where there 
is only one nearby galaxy, and no other obser- 
vational constrains, such as a second image, we 
assume a Singular Isothermal Sphere (SIS) mass 
profile for the lensing galaxy. This model predicts 
multiple images where the impact parameter, 0^, 
is smaller than twice a galaxy's Einstein radius, 
6e, which is defined as 



= 47r 



(^) 



2 A. 



(1) 



Here ai, is the line-of-sight velocity dispersion of 
the galaxy (which could be constrained from the 
TF analysis described above), and Dis and Ds are 
the angular diameter distances between the lens 
and the source, and the observer and the source, 



respectively (Schneider et al. 1999). In addition 



to lensing geometry, the model also provides in- 
formation regarding intensity magnification and 
time delay between images. In the simplest case 
of a single lensing galaxy, the first image is further 
from the center-of-mass of the lens and is brighter 
than the second image, making it more likely to be 
detected. More complex lensing systems are not 
bound by such constraints, but still make predic- 
tions about magnifications and time delays. 

Although Mgll absorbers and intervening 
galaxies are correlated for quasar sight-lines, 
the relation between the impact parameter and 
the Mgll equivalent width is not a tight one 



(Churchill et al. 2005), making the task of iden- 



tifying the right absorber for each system not 
straight-forward. As we normally lack the data to 
verify the nearby galaxies' redshifts, we implement 
a statistical test to understand the likelihood of 
a nearby galaxy being at the absorber's redshift. 
Again using the UDF galaxy catalog of |Coe et al.| 
(2006), we perform a Monte-Carlo simulation that 



measures the probability of having a galaxy with 
a given apparent magnitude within a known dis- 
tance from a random line of sight. This allows us 



to estimate the chances of finding galaxies near a 
GRB. 

Hereafter, we analyze the possibility of strong 
gravitational lensing for each of our GRBs, in de- 
scending order of multiple-lensing likelihood. 

4.1. Candidates of Multiply Imaged GRBs 

4.1.1. GRB020405 

HST images of this GRB (z=0.695) revealed 
another transient object 3" away, as noted by 



Masetti et al. (2003). A VLT spectrum confirmed 



that both the nearest galaxy and one other (ob- 
jects 1 and 2, respectively, in Figure [2| are at the 
redshift of the z=0.472 Mgll absorbing system. 
The second transient is visible along the edge of 
galaxy 1. Objects 3,4,5 and 6 are nearby galaxies 
which we speculate below could possibly be part 
of a group at the absorber redshift. 

The second transient was observed in the first 
observation of HST, 23 days after the GRB trig- 
ger, and was fainter than the GRB at all times. 
With only few data points for this transient, it is 
consistent (within errors) with the expected colour 
evolution of the GRB. While we cannot dismiss 
the possibility of the second transient being a non- 
related field SN, we expect the probability of find- 
ing an object 3// away from any z ^ 0.5 redshift 
galaxy is low. Modelling this system with LENS- 
MODEL, we are able to explain the unknown tran- 
sient as an earlier, brighter image of the GRB 
(Rapoport et al. in prep). The model assumes 
galaxies 1-6 are a group at the Mgll absorption 
redshift. We assume the galaxies lie in common 
dark matter envelop and use a Singular Isother- 
mal Ellipse (SIE) mass profile with shear. While 
the time delays, positions of images and their mag- 
nification are highly sensitive to the model, pos- 
sible solutions include a time difference between 
the images of ^120 days, with the first GRB im- 
age being 1.8 times brighter. One such model is 
shown in figure [5] The transient (i.e., first GRB 
image) would not have been detected in optical 
at the time of the original GRB020405 observa- 
tions, as the resolution and depth required to sep- 
arate the transient from the nearby galaxy were 
not achieved by the ground-based facilities used 
in the first few weeks after the gamma-ray trigger. 
Figure [3] shows the GRB observations in I band in 
black and the predicted fiux for the first image in 
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red (1.8 times brighter). The blue points are the 
observed fluxes of the second transient from the 
HST images and the upper limit was found by sub- 
tracting the last HST observation from the VLT I 
band image. (The earlier epochs for GRB020405 
were obtained with smaller telescopes, and pro- 
vide no additional constraints on the light curve 
of the second transient.). We are pursuing further 
observations to verify the redshift of the galax- 
ies in the field. The model predicts several less- 
magnified images of the host complex, which do 
not conflict with current observations. In addi- 
tion, some models suggest there were even earlier 
images of the GRB (on the right of Fig. [2]), which 
occurred years before. It is not possible to explain 
the second transient as a GRB image if the nearby 
galaxies are not at the absorbers redshift. Phot- 
Z analysis, using the EASY software (Brammer 
et al.||2008[ ), of the possible group members shows 



consistency for the galaxy redshifts, but with large 
uncertainty. 

4.1.2. GRB030429 

The impact parameter between the GRB030429 
and the nearby galaxy is 1.2" (flgure|4]), and the 
redshift of the galaxy was confirmed by [Jakobs 



son et al. ( 2004 ) to be that of the absorbing Mgll 



system. Moreover, at a redshift of 2.66 with an ab- 
sorbing system at a redshift of 0.8418, the angular 
diameter distance ratio is approximately 0.54. 
As gravitational lensing is most efficient when the 
angular size distance of the source as seen by the 
lens equals the angular size distance between the 
lens and the observer, GRB030429 represents a 
likely candidate for being lensed. 

Our Monte-Carlo simulation with the UDF 
finds the probability of having a galaxy with the 
observed magnitude within 1.2" along a random 
line of sight to be 0.4%. For the lines of sight 
towards our 11 GRBs, the chances of randomly 
finding one such alignment is ^14%. Using the 
magnitudes of the galaxy given by |Jakobsso"n| 
et al. (2004), our SED fitting techniques imply 
Mb = -21.1 ± 0.1. The TF relation for such 
an intrinsically bright galaxy suggests a velocity 
dispersion of 160 ± 65 km/s which corresponds to 
9e = 0.38to;25 arcseconds. With a separation of 
1.2 arcsec, the SIS model predicts a second image 
for this GRB only if <; 200 km/s. Since lensing 
would select out those galaxies with the highest 
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Fig. 2.— Left: HST WFPC2/F702W field of 
GRB020405. The GRB is clearly visible and the 
complex host and second transient are indicated. 
Objects 1 and 2 are confirmed to be at the strong 
absorber's redshift of z=0.472. Galaxies 3-6 may 
be part of a group at the absorber's redshift. 
Right: A representative LENSMODEL solution 
for SIE with shear. The time delay between the 
GRB images is ~120 days where the leading im- 
age is ~1.8 times brighter that the later one. The 
two predicted GRB images on the right arrive hun- 
dreds of days earlier. The host images to the right 
are less magnified. 



(Tj/, the TF relationship is useful to show plausi- 
bility in this case. Because the Einstein radius 
for this system needs to be ^0.6" for multiple 
imaging, the impact parameter of the observed 
GRB means it was likely the first of two images 
if lensed. Figure [5] shows the expected magnifica- 
tions of the GRB and the predicted second image, 
and the time delay between the two. The last 
observation of this GRB was taken with the VLT 
~67 days after the trigger, which would have been 
too soon for detecting the second image, assuming 
a^, < 250 km/s. No late X-ray observations were 
taken. In order to test this being a lensing system, 
direct measurements of the velocity dispersion or 
galaxy mass are needed. 
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Fig. 3. — Light curve prediction for the lensing 
model. In black dots, the observed I band flux of 
GRB020405. In red, the predicted light curve of 
the first image, magnified by 1.8 with a time delay 
of -120 days. The blue dots are the observed flux 
of the second transient in the HST images, and the 
blue upper limit is found when subtracting the last 
HST image from the early VLT I image. 




Fig. 4.— VLT/F0RS2 R_BESS field of 
GRB030429, with the potential lensing galaxy 
-1.2" to the right. 



4.2. Possible Candidates of Multiply Im- 
aged GRBs 

4.2.1. GRB021004 

While there is no exceptionally bright galaxy 
near this GRB, there are five faint galaxies within 
4.5" which could be a part of a group (Fig.|6]). At a 
redshift of 2.33, this GRB had 2 strong absorbing 
systems at redshifts 1.38 and 1.60. 
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Fig. 5. — The prediction from a SIS model for 
the magnification of GRB030429 (blue), its second 
image (green), and the time delay between them 
(red). 



Fynbo et al. (2005) studied the colors of the 
neighboring galaxies and found the host galaxy 
to be different from its surrounding. As gravita- 
tional lensing does not alter the color of objects, 
this suggests that none of the other galaxies is a 
second image of the host. Therefore, if modeling 
the group as an SIE and for the case of a dou- 
ble image, the detected GRB must have been the 
first, brighter one, and the second image of the 
host galaxy must not be strongly magnified. Us- 
ing the GRAVLENS software we construct a pos- 
sible model of the group, having a common dark 
matter envelope which can be approximated as a 
SIE. The model predicts multiple images, with the 
second being too faint to detect the host. The cen- 
ter of mass of the group model appears towards 
the more dense part of the group, as expected. 
Since the time delay between the images is ~550 
days, we would not expect to observe 2 images of 
the GRB at the same time. Unfortunately, the 
HST observations stopped 53 days after the trig- 
ger, which is less than the expected delay time for 
this model. The host galaxy observations taken 
a year later would have still been to early. Al- 
though this GRB had a bright optical afterglow, 
with a relative demagnification of a factor of —25, 
the window of possible detection is fairly narrow 
(several days). 

Redshift confirmation for the galaxies is re- 
quired to verify this is indeed a bona fide group. 
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In such a case, deeper images might reveal other 
background sources and possibly multiple images 
of the host. Phot-z analysis using the HST im- 
ages finds galaxies 1-5 to be consistent with the 
absorber's redshift within errors. 
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Fig. 6.— Top: HST ACS/F606W field of 
GRB021004. The host complex and GRB are 
marked, and object 1-5 are the neighbouring 
galaxies. Bottom: Possible GRAVLENS solution 
for GRB021004. The galaxy group was modeled 
as a SIE. The observed GRB is ^25 times brighter 
than the second image, and arrives ^550 days ear- 
lier. 



4.2.2. GRB010222 

A strong Mgll absorbing system was found at 
a redshift of 0.927 for this z=1.477 GRB. While 
there was no single bright, close galaxy, HST im- 
ages show a crowded region with at least 11 clearly 
identified galaxies which could, again, be a part of 
a group (Fig. [t]). To produce a second image, the 
velocity dispersion of the group would need to be 
250 km/s for a group centered 0.5" from the GRB, 
350 km/s for 1" between the GRB and group cen- 
ter, or 500 km/s for 2". 
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Fig. 7.— HST/WFPG2 F606W image of 
GRB010222. The host galaxy of the GRB is 
clearly visible. At least 11 galaxies are seen in 
the image, and could possibly be a part of a group 
(objects 1-11). 



Here again, one would need to confirm the 
galaxies are at the absorber redshift as the first 
step to concluding that this group lensed the GRB. 
As the SNR of these galaxies is low and the dumpi- 
ness of some making it difficult to measure accu- 
rate photometry, we performs a phot-z analysis 
on galaxies 1-4,9-11. The results, whilst carrying 
large errors are consistent with the galaxies being 
at the absorber's redshift. 
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4-2.3. GRB060206 

At a redshift of 4.048, GRB060206 exhibits 2 
strong absorbing systems at redshifts 2.26 and 
1.48. An HST image reveals 2 galaxies near the ob- 
served host (see Fig. [s]), at distances of ~0.96" for 
galaxy 1 and ^ 2.45" for galaxy 2. The probability 
of finding an alignment of two such galaxies along 
a random line of sight is 3%, or a 29% chance of be- 
ing observed in the 11 GRBs. The field was also 
observed with CMOS r' on Gemini North (GN- 
2007A-Q-88), and galaxy 2 was detected with a 
magnitude 23.9±0.1 (|Thone et al]|2008[). Using 



the color information for galaxy 2, we fit an SED 
model and estimate it could be a small starburst 
galaxy at either absorber redshift, with Mb ^ - 
16.35 (-17.46) for z=1.48 (2.26). Galaxy 2 would 
have to have a velocity dispersion larger than 325 
km/s to create 2 images of this GRB. Since we do 
not expect such velocities from a small starburst 
galaxy, we do not consider the galaxy further. 

With galaxy 1 only observed with one filter, it is 
impossible to determine the galaxy type. Galaxy 1 
would require an Einstein radius of ^ 0.5" to lens 
the GRB, which is equivalent to cr ~ 200 (275) 
km/s at redshift 1.48 (2.26). A second image of 
the GRB host would not necessarily be observed 
due to the small separation from the lensing galaxy 
and the usual de-magnification of the second im- 
age relative to the first image. Further photom- 
etry is required for constraining the model of the 
galaxy and determining its lensing feasibility. 




Fig. 8.— HST ACS/F814W field of GRB060206. 
Object 1 (2) is - 0.96" (2.45") from the host 
galaxy. 



4.3. Unlikely Candidates of Multiply Im- 
aged GRBs 

4.3.1. GRB060418 

This GRB at z=1.49, with 3 strong absorbing 
systems at z=0.603, 0.656, and 1.107, was well 



studied by Pollack et al. (2009), who identified 



a complicated host and 3 adjacent galaxies (fig- 
ure |9]). In their paper, they identify galaxy 1 
as that responsible for the absorbing system at 
z=0.656, as its redshift was confirmed with ob- 
served emission lines. Object 2 was assumed to be 
the z=1.107 absorber and object 3 the z=0.603 ab- 
sorber due to their angular sizes. However, none 
of these later identifications were confirmed via 



spectroscopy. Recently, Chen (20111 studied the 



nearby galaxies and was able to spectroscoply con- 
firm the redshift of what is now identified as the 
host complex in figure |9] 

Due to the relatively large distance between the 
GRB and the probable absorbers, we find lensing 
to be unlikely. 




Fig. 9.— HST/WFPC2 F775W image of 
GRB060418. The GRB and host locations are 
marked, 1-3 are nearby galaxies. Galaxy 1 was 
confirmed at the absorber z=0.656 redshift. 



4.3.2. GRB050820A 



In their paper, Chen et al 
GRB at z 
at z 



(20091 study this 



2.615 and its strong absorbing systems 
0.6915 and 1.4288. They identify 2 compact 
objects at separations of 1.3" and 0.4" from the 



GRB. In spectra taken lately by Chen (2011 1, the 



objects were found to be a part of the GRB host. 
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Therefore, we conclude that this GRB was not 
strongly lensed. 

4.3.3. GRB080319B 

This GRB at a redshift of 0.9378, with one 
strong Mgll intervening system at z=0.7154 and 3 
other weak systems, has 3 nearby galaxies. Find- 
ing 3 galaxies within 4" is unlikely, with the mean 
number of galaxies predicted to be 1.2, with a 
standard deviation of 0.45. However, the faint- 
ness of the galaxies makes it impossible to find a 
model which allows lensing. Although these data 
were taken with only 2 filters, we were able to fit 
an SED model which found them to be faint early 
type galaxies, with —15.5 < Mb < —12.5 (as- 
suming z=0.7154). Therefore, even if the galax- 
ies are in the absorber's redshift and not back- 
ground galaxies, they could not be massive enough 
to cause strong lensing. 

4.3.4. GRB991216 




Fig. 10.— HST/STIS Long Pass image of 
GRB991216. The GRB host and the nearby ir- 
regular galaxy are marked. 

A faint (R^24.5) irregular galaxjj^lies 0.4" to 
the South of this GRB. The GRB redshift was 
1.022 and the intervening systems were at 0.77 
and 0.803. This configuration requires that nearby 
galaxy to have a velocity dispersion of at least 190 
km/s, which is unlikely for a relatively faint irreg- 
ular galaxy. However, the proximity implies that 
even a velocity dispersion of 70km/s would imply 
15% magnification. Therefore, while not doubly 

5 http://www.stsci.edu/~fruchtcr/GRB/991216/index.html 



lensed, this GRB (and its host) are likely to be 
magnified if the nearby galaxy is at the absorber's 
redshift. 

4.3.5. GRB020813 

At a redshift of 1.255 this GRB shows a strong 
MgH absorber at a redshift of 1.224. The ratio 
of is ^^58 which suggests that strong lensing 
is highly unlikely. The nearby bright galaxy at 
a distance of 2.4" away would require > 1550 
km/s, which is unphysical. 

4.3.6. GRB050908 

This GRB at z=3.55 with a strong MgH 
absorbing system at z= 1.548, reveals a faint 
(F775W~'26.6 mag) galaxy ^ 1" away. This align- 
ment would require the galaxy to have a velocity 
dispersion larger than 230 km/s. HST observa- 
tions with different filters are required to model 
the SED of this galaxy and conclude if this high 
velocity dispersion is feasible. 

5. Discussion and Conclusions 

There is a well known excess of strong MgH ab- 
sorbers towards GRBs compared to quasar lines 
of sight. The most viable explanations for this 
discrepancy are gravitational lensing of the GRB 
population, and dust extinction of the population 
of quasars with strong absorbers. In this paper 
we have identified a new difference between GRBs 
and random lines of sight. We find that GRBs 
with strong MgH absorbers are found closer to 
bright galaxies when compared with random lines 
of sight that are restricted to the 60% of sky near- 
est to foreground galaxies (corresponding to the 
fraction of GRBs with strong MgH absorption). 
This new property of the GRB population can- 
not be explained by quasar extinction. On the 
other hand such a concentration of GRBs around 
foreground galaxies is a natural consequence of 
gravitational lensing. We therefore argue that 
strong gravitational lensing is likely to be play- 
ing a role in explaining the discrepancy between 
MgH absorbers towards GRBs and quasars. |Por-| 
ciani et al. also found that GRBs with MgH ab- 
sorbing systems were slightly brighter (a factor of 
1.7) than those without, implying again that grav- 
itational lensing could explain the phenomenon. 
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provided the GRB luminosity function was suffi- 
ciently steep. 

Wll predicted that 10%-60% of the 26 GRBs in 
the V09 sample should have been multiply imaged 
if lensing is the explanation for the enhanced Mgll 
absorption towards GRBs relative to quasars. As- 
suming that the GRBs with strong Mgll systems 
(EW> lA) are the only ones which are multiply- 
imaged, we would expect 1.9 - 11 of the 11 GRBs 
studied to be strongly gravitationally lensed. If 
none of our sample were plausibly gravitationally 
lensed, then this would effectively rule out the 
gravitational lensing hypothesis. 

Quasars are observed to show strong lensing 
~0.1% of the time, and we expect, if the GRB 
and quasar samples are similarly distributed, no 
GRBs in our sample to be lensed. It is important 
to note that the Wll model predicts a lensing rate 
of less than 1% for both GRBs selected only in 
gamma-rays and for quasars. Thus, even the null 
detection of multiple imaging in the 2700 GRBs 
observed by the BATSE instrument, which could 
see a much larger fraction of the sky than Swift at 



any given time, is not surprising (see Porciani & 
Madau||2001|). 



Our analysis reveals a high frequency of galax- 
ies near GRB sight-lines, with 2 suspected two- 
image systems (GRB020405 and GRB030429), 2 
possible galaxies groups that would result in mul- 
tiple images (GRB021004 and GRB010222), and 
several other less probable candidates. Since only 
the two most likely candidates have spectroscopic 
confirmation that the nearby galaxy lies at the ab- 
sorber's redshift, more data are needed to confirm 
the strong gravitational lensing hypothesis. This 
study, while not conclusively demonstrating lens- 
ing, is consistent with what would be expected 
from the model of Wll. If, as described by Wyithe 
et al., the cumulative luminosity function of the 
optical afterglow follows a power law of the form 
\E'a(-^^a) oc (Lyi)""-*, where La is the afterglow 
luminosity, then our study indicates a value of 
aA ~ 3.5 or steeper. 



Recently, Budzynski & Hewett (2011) studied 



the correlation between quasars, Mgll absorbers 
and extinction, and found strong dependence of 
E(B— V) on the absorber's equivalent width. Mod- 
eling the observed difference between GRBs and 
quasars, they find dust obscuration to be a signifi- 
cant factor in the Mgll discrepancy at high equiv- 



alent width. However, they acknowledge that it 
is unlikely to be the only effect, and that the full 
explanation involves more than one process. Our 
analysis compares the galaxy over-density towards 
GRBs relative to a random line of sight and so 
is independent from comparisons to other objects 
(e.g., quasars or blazars; see Bergeron et al.|20lT |. 
Therefore, both lensing and dust could be of sig- 
nificance for solving the Mgll problem. 

Searching for galaxies near the lines of sight to- 
ward quasars with strong Mgll absorption would 
require detailed modeling and subtraction of the 
quasar point spread function, which is beyond the 
scope of this paper. The virtue of the GRB follow- 
up is the fading of the point source, which affords 
a straight-forward assessment of the surrounding 
sky. If follow-up studies of the most probable 
lenses do not unambiguously rule in or out lens- 
ing, then a future direction could include apply- 
ing a similar analysis to a quasars sample. How- 
ever, in order to best test the gravitational lensing 
model, what is needed is a real time study. Once 
an appropriate candidate is identified, deep high- 
resolution images should follow, to search for the 
possible lensing system. If one is found, lensing 
models should be applied to predict if a second 
image is expected. If so, the location and time de- 
lay of the next image can be measured. Depending 
on the error expected in the time delay, which is 
proportional to crjj, an appropriate observational 
cadence can be put into action. Since the sec- 
ond image is the inner one for the simple SIS/SIE 
models, it is likely to suffer from stronger extinc- 
tion by the lensing galaxy. Moreover, the rapid 
fading in optical and the smaller magnification of 
the second image compared to the first one might 
require frequent observations in order to catch the 
second image while it is observable. Therefore, we 
suggest conducting the search in X-rays, which are 
less affected by passing through a galaxy. 

The preferential lensing of GRBs would herald 
a shift in the study of both gravitational lenses and 
GRBs. Prior knowledge of where and when a GRB 
image will appear should allow coordinated obser- 
vations from earth and space to be scheduled, lead- 
ing to unprecedented multi- wavelength data. Such 
information would be of great assistance in under- 
standing progenitor and outflow properties. More- 
over, accurate measurements of the time delays 
between lensed images would constrain cosmolog- 
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ical models and inform studies of the dark matter 
distribution associated with the lensing galaxy. 

We hope that the conclusions drawn in this 
work will stimulate the search for lensed GRBs 
(perhaps identified by their strong Mgll ab- 
sorbers), and the associated tertiary images. 
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